Inefficient emergent oscillations in intersecting driven many-particle flows 
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Oscillatory flow patterns have been observed in many different driven many-particle systems. The 
conventional assumption is that the reason for emergent oscillations in opposing flows is an increased 
efficiency (throughput). In this contribution, however, we will study intersecting pedestrian and 
vehicle flows as an example for inefficient emergent oscillations. In the coupled vehicle-pedestrian 
delay problem, oscillating pedestrian and vehicle ffows form when pedestrians cross the street with 
a small time gap to approaching cars, while both pedestrians and vehicles benefit, when they keep 
some overcritical time gap. That is, when the safety time gap of pedestrians is increased, the average 
delay time of pedestrians decreases and the vehicle fiow goes up. This may be interpreted as a slower- 
is-faster effect. The underlying mechanism of this effect is explained in detail. 

PACS numbers: 89.40.+k,47.54.+r,89.75.Kd,47.62.-|-q 



Emergent oscillations have been discovered in so dif- 
ferent systems as the density oscillator 0, ticking hour 
glass 1^, RNA Polymerase traffic on DNA Isl, pedestri- 
ans passing a bottleneck |J, l5| , or ants 6] . Despite the 
different underlying oscillation mechanisms and system 
components, all of these systems can be treated as driven 
or self-driven many-particle systems characterized by 
counterflows. Therefore, one might assume a unified prin- 
ciple behind these emergent oscillations such as an opti- 
mization of throughput: The clustering of units with the 
same flow direction could reduce "frictional" interactions, 
which are particularly high between units moving in dif- 
ferent directions. In this contribution, we will study the 
example of intersecting pedestrian and vehicle flows. This 
system is found to show a transition to emergent oscil- 
lations as well. However, contrary to our expectations, 
the oscillations are not an efficient pattern of motion. In- 
stead, they are related with a considerable reduction of 
the throughput and increased waiting times. 

In the past, the investigation of vehicle and pedestrian 
streams by means of experiments and models from sta- 
tistical physics or fluid-dynamics has revealed the mech- 
anisms behind many observedphenomena such as differ- 
ent forms of congestion [dJMJSj. Moreover, various self- 
organization phenomena [J, 1^ have been discovered in 
pedestrian flows, including the so-called "faster-is-slower 
effect" in "panicking" crowds j^. These have stimulated 
research in many other fields such as colloidal 0| and 
biological systems. 

Let us now come back to the problem of interacting ve- 
hicle and pedestrian flows, a problem that has not been 
thoroughly studied in the past. In a way, the problem 
can be viewed as two dynamically coupled queues, which 
cannot be served simultaneously, since pedestrians must 
cross the street at times when no vehicle passes and vice 
versa. Such coupled queuing systems are known to dis- 
play interesting dynamic behaviors, including irregular 



oscillations jl^] and chaos [ijl- We are, therefore, inter- 
ested in identifying the possible dynamic behaviors of 
the coupled vehicle-pedestrian-delay problem, their per- 
formance and preconditions. 

The pedestrian delay problem is a growing concern of 
urban planning. It is defined as follows: Suppose there 
is a stream of vehicle traffic moving on a main street, 
and suppose that a pedestrian arrives at time t — to at 
the roadside and intends to cross this street (away from 
any pedestrian crossing facility), see Fig. 1(a), (b). What 
is then the average delay to the pedestrian? The early 
pedestrian delay models assume that there is a negative 
exponential distribution of vehicular headways 01 EB| ■ 
Other models have adopted a shifted exponential distri- 
bution, a double-displaced negative exponential distribu- 
tion, etc.pjj. Recently Guo et al. yjl have proposed a 
pedestrian delay model, in which the overall delay to 
pedestrians is obtained as a combination of the delay 
by traffic-light induced vehicle clusters and the delay to 
pedestrians arriving during the random vehicle flow be- 
tween the clusters. 

Note that, in the pedestrian delay problem, the interac- 
tions between vehicles and pedestrians, at least influences 
of pedestrians on vehicles, have not been considered, yet. 
Usually, it is assumed that the crossing of pedestrians will 
not affect the motion of vehicles. This is certainly not re- 
alistic. Therefore, this paper studies the coupled vehicle- 
pedestrian delay problem, taking into account mutual 
interactions. This is relevant for the capacity of traffic 
infrastructures for both, vehicles and pedestrians. 

Our delay model is as follows: Firstly, within one in- 
cremental time step dt of 0.1s (corresponding to the ap- 
plied time discretization of the car- following model), we 
assume the arrival of one pedestrian along the roadside 
with probability p. When a pedestrian arrives along the 
roadside at a given point O, he or she checks the traffic 
situation [Fig. 1(a)]. We will distinguish two situations: 



2 



(i) No other pedestrian is on the road. In this case we 
suppose that, when the safety criterion 



d > do + at-cVn 



(1) 



is satisfied, the pedestrian will cross the road. Here, do is 
the minimum safety distance of pedestrians, the time 
needed for a pedestrian to traverse a one-lane street, tr a 
safety coefficient, d the distance from the nearest vehicle 
n upstream of point O, and Vn its velocity [13. (u) Other 
pedestrians are crossing the road. In this case, the pedes- 
trians on the road will encourage newly arriving pedestri- 
ans to follow, as an obstructed driver-vehicle unit would 
not dare to accelerate. This effect can be simulated by 
adjusting the safety coefficient a. 

We assume that in case (i), the safety coefficient a 
chosen by a pedestrian is (Tq: while in case (ii), he or she 
will choose a smaller safety coefficient cti . The pedestrian 
simulation is similar to that described in Ref. J19|: First, 
the position Xn{t) of the nearest vehicle n upstream of 
the crossing point O is identified. If a pedestrian is on the 
street, the net distance to the next object is specified as 
Ax{t) = d{t) = xo — Xn{t), and the velocity of the object 
ahead into the driving direction is w = 0. Otherwise the 
distance and velocity are given by the next vehicle n — 1 
ahead, i.e. Ax{t) — Xn-i{t) — Xn — and v{t) = 

Vn-iit)- This enters the equation of vehicle motion 



dVn 

~dr 



f{Ax,V,Vn) + C„(i) . 



(2) 



Here, / is the vehicle length, / the acceleration function 
and ^ a stochastic term. For illustrative purposes, the 
acceleration function has been specified according to the 
well-investigated intelligent driver model (IDM) 



f{Ax,v,Vn) = a 



Aa; 



(3) 



The parameter vq denotes the desired velocity, while s* 



So -I- TVr, 



is the desired minimum gap, where 



So is the minimum safety distance of cars, T is the safe 
time gap, a the maximum acceleration, and b the desired 
deceleration. The stochastic term ^ has been set to zero. 

In each simulation time step, if a random number y 
uniformly distributed between and 1 is smaller than 
p, a new pedestrian occurs at point O. The pedestrian 
crosses the road if the safety criterion |^ is satisfied. 
Otherwise, he or she will wait until next time step. 

In the simulations, we adopt open boundary condi- 
tions. Suppose that a car has been added at time t = Iq, 
then another car is added at place x = and time 
t = to -\- At, provided that the previous car is at least 
7 m away (2]|. Here At is determined by a shifted ex- 
ponential distribution: At — ^^-[jt^t + 2 s, where z is 
a uniformly distributed random number between and 
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FIG. 1: (a), (b) Sketch of the vehicle-pedestrian delay system, 
(c) Average delay time of pedestrians and (d) average vehicle 
flow for (Ti — 0.6, both as a function of the pedestrian arrival 
probability p. 



1. For simplicity, the speed of the new car is assumed to 
agree with the one of the car ahead. At the exit {x = L), 
cars are removed. If there is no leading car, the distance 
Ax is set to some large number and v to the maximum 
velocity vo- 

In the simulations, Eq. jSJ is solved by the Euler 
method. The time step is set to 0.1 s, as smaller val- 
ues do not change our results. The model parameters are 
Wo = 15 m/s, a = 2 m/s^, b = 1.5 m/s^, T = 1.4 s, 
So = 2 m, = 2 s, Z = 5 m, do — 1.2 m, and ao = 2. 
Note that a probably less realistic model with fewer pa- 
rameters could be used as well. 

We first show simulation results for cti = 0.6. In our 
simulations, vehicles enter the empty road from t — on 
and the first pedestrian arrives after t = 500 s. The road 
length is L — 1400 m and the location of the crossing 
point xo = 1200 m. 

Figure 1(c) shows the average pedestrian delay. One 
can see that, when the arrival probability p is small, the 
average delay time essentially remains constant. How- 
ever, when p > 0.05, it begins to increase with grow- 
ing values of p. Then, after reaching the maximum at 
p w 0.32, it goes down with a further increase of the 
arrival rate r = p/dt. 

Next, let us try to explain the change of the slope. For 
this, we show in Fig. 2(a) the typical structure of the 
traffic situation obtained aX p — 0.25. One can see that 
at the point O {x — 1200 m), there is an alternating 
vehicle and pedestrian flow. When pedestrians cross the 
road, the vehicles are stopped. When the stopped vehi- 
cle queue starts [this occurs when no pedestrian arrives 
in the subsequent 20 time steps (corresponding to the 
crossing time ij- = 2 s)], the formed vehicle queue does 
not allow pedestrians to cross until a large gap occurs, 
again. 

Figure 3(a) sketches the underlying mechanism. If 
pedestrians have stopped a vehicle at time ^O: the fol- 



3 




1000 1200 1400 



1000 1200 1400 



X (m) X (m) 

FIG. 2: (a) Representative space-over-time plot of vehicle tra- 
jectories for (Ji = 0.6. Sometimes, vehicles are stopped by 
crossing pedestrians and form a queue. This suppresses the 
crossing of newly arriving pedestrians for a long time period 
later on. (b) Representative space-over-time plot of vehicle 
trajectories for cri = 1.0. Obviously, vehicles are decelerated, 
but not stopped by pedestrians. 
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FIG. 3: (a) Sketch of the queue formation and queue res- 
olution mechanism (see text), neglecting details of acceler- 
ation and deceleration, (b) Dynamically changing time gap 
[d(t) — d(^lv{€) of a vehicle if pedestrians enter the street 100 
meters ahead (see text). If crifr < 2, pedestrians continue en- 
tering the street, which may stop the vehicle (solid growing 
curve). Otherwise, the crossing criterion is violated after some 
time and the vehicle can accelerate (dashed falling curve). 



lowing vehicles queue up with a speed of |2S 



C 



1 



(4) 



and more pedestrians can cross, pjam is the jam density, 
Qarr the arrival rate of vehicles and uq their desired ve- 
locity. After a time interval A^i, i.e. 20 time steps after 
the last pedestrian has entered the road, the first vehicle 
in the queue can accelerate again. The last vehicle in the 
queue reaches point O at time to + + with [i^l 



At2 = CAii 



1 + \c\lvQ 
\c\-C 



(5) 



where c ~ "l/(pjamr) w —15 km/h is the characteristic 
jam resolution speed. Afterwards, pedestrians have the 
chance to find a gap in the vehicle flow, again. 



During the period of , the average number of arriv- 
ing pedestrians is p At2. We can assume that the waiting 
time of the last pedestrian is approximately zero, while it 
is approximately Ai2 for the first one. The average delay 
is At2/2, and the overall expected delay of all waiting 
pedestrians amounts to p(At2)^/2. As the delay to the 
pedestrians arriving during the time period Ati is zero, 
the average delay is 



p(At2)V2 + 0-pAti 
p(Ati + A<2) 



(At2)^ 



2(Ati + At2) 



(6) 



With the increase of p, the probability (1 — that 
no pedestrian arrives in the 20 subsequent time steps 
decreases. This implies that the average value of Aii goes 
up, which leads to an increase in the average delay. 

Nevertheless, when p is large, the stopped vehicle 
queue may occupy the whole length of the road upstream 
of point O. In this case, At2 will become a constant. 
Therefore, with a further increase of p, the average delay 



(At2)^ 



decreases, as Ati increases with p |23 



2(Ati-|-At2) 

Let us now have a look at the average traffic flow 
of vehicles shown in Fig. 1(d). One can see that, when 
■p < 0.28, the vehicle flow remains essentially constant. 
When p > 0.28, it decreases with an increase of p. This is 
because a stopped vehicle queue may occupy the whole 
street upstream of point O for a certain time period, so 
that the entry of further vehicles becomes temporarily 
impossible. In other words, the crossing pedestrian flow 
affects the street capacity for p > 0.28. 

We will now increase the safety coefficient of pedestri- 
ans to tTi = 1.0. This introduces an additional safety time 
gap for crossing the road and, thereby, reduces the im- 
pact on approaching vehicles [see Fig. 2(b) for p = 0.25]. 
Figure 1(c) shows simulation results for ai = 0.6 and 
CTi = 1.0. One can see that, compared with the result 
for CTi = 0.6, the average delay time is slightly higher at 
small values of p, but it is considerably smaller for large 
arrival probabilities p. 

In order to understand this qualitatively different re- 
sult, we have studied the deceleration process of a freely 
moving vehicle when a driver sees pedestrians entering 
the street Aa; = d = 100 m ahead, with arrival probabil- 
ity p = 1. In Fig. 3(b), the solid line shows the evolution 
of the time gap [d{t) — dQ\/v{t) with time t in the case of 
CTi — 0.6. First, the time gap decreases with time. How- 
ever, after it reaches a minimum value, it increases again. 
During the whole deceleration process, the time gap is 
larger than aitr- This means that pedestrians can always 
enter the road, and the vehicle may be stopped by pedes- 
trians. However, for cti ~ 1.0, pedestrians will not enter 
the road after the time gap has decreased to tritr = 2 s 
at time to, as the safety criterion becomes violated. 
Nevertheless, the vehicle continues decelerating during 
the crossing time t,- = 2 s. Afterwards, the vehicle ac- 
celerates again [dashed falling curve in Fig. 3(b)] and no 
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pedestrian can enter before the vehicle has passed point 
O. This means that pedestrians will never stop vehicles. 
Consequently, no vehicle queue will form and pedestrians 
will cross the street one by one or in small groups [see 
Fig. 2(b)]. Our simulations show that the transition be- 
tween the continuous and oscillatory crossing dynamics 
occurs at cti w 0.96. 

We have investigated the coupled pedestrian-vehicle 
delay problem with a simplified model. In contrast to 
our approach, previous pedestrian delay models have ne- 
glected the influence of pedestrian crossing on vehicle dy- 
namics. Our computer simulations were carried out for 
two different types of pedestrians to highlight the transi- 
tion of the system behavior at ai « 0.96: (i) Aggressive 
pedestrians with a small safety coefficient cti < 0.96 could 
force vehicles to stop by successive crossing events, which 
produced alternating vehicle and pedestrian flows at high 
pedestrian arrival rates r = p/dt. (ii) Careful pedestrians 
with a safety coefficient cti > 0.96 did not stop vehicles 
and crossed the street one by one or in small groups. 
Altogether, this mode was more efficient for pedestrians 
and cars, as single pedestrian could not keep a growing 
vehicle queue from going. 

As in panicking crowds 0, impatience affects system 
performance in a negative way: Waiting a bit longer 
(for a larger vehicle gap) implies smaller average de- 
lays ( "slower-is- faster effect"). Therefore, when the ve- 
hicle flow is not too large, a traffic light is not needed 
to allow pedestrians to cross the street. It is rather re- 
quired to terminate inefficient pedestrian crossing while 
vehicle queues are building up |24| . Future investigations 
should clarify, whether the emergent oscillations in the 
density oscillator, ticking hour glass, pedestrians streams, 
ant traffic, and collective motion of molecular motors are 
also inefficient, as in the surprising example discussed 
here, or efficient, as expected. 
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